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The systematic theory of the formation of the short light pulses in the squeezed state during the 
propagation in a medium with inertial Kerr nonlinearity is developed. The algebra of time-dependent 
Bose-operators is elaborated and the normal-ordering theorem for them is formulated. It is estab- 
lished that the spectral region where the quadrature fluctuations are weaker than the shot-noise, 
depends on both the relaxation time of the nonlinearity and the magnitude of the nonlinear phase 
shift. It is also shown that the frequency at which suppression of the fluctuation is greatest can 
be controlled by adjusting the phase of the initial coherent light pulse. The spectral correlation 
function of photons is introduced and photon antibunching is found. 

PACS numbers: 42-50, 42-50. L, 42.50.Dv 
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I. INTRODUCTION 

During the past years the study of formation of non- 
classical light pulses in nonlinear media has been the fo- 
cus of a considerable attention. The present article is 
devoted to the development of the consecutive theory 
of formation of nonclassical short light pulses in non- 
linear media with inertial Kerr nonlinearity. It is well 
known, that in nonlinear media in the presence of the 
self-action effect the squeezing of quantum fluctuations 
of one quadrature component of a field with conserva- 
tion of the photon statistics take place 0. At present 
there are two basic directions of research in the quantum 
theory of self-action of ultrashort light pulses (USPs). 
In the first approach [^-|| the calculations of the non- 
classical light formation at the self-action of light pulses 
assume that the nonlinear response of nonlinearity is in- 
stantaneous and that the relative fluctuations are small. 
This latter assumption is valid for the intensive USPs 
frequently used in experiments. However, a finite relax- 
ation time of the nonlinearity has a principled impor- 
tance as the relaxation time determines the region of the 
spectrum of the quantum fluctuation below the standard 
noise level. For the first time, in || was noted that for 
the correct quantum solution of self-action, it is neces- 
sary to take into account the presence of quantum noise. 
The presence of quantum noise was anticipated in |?j as 
thermal addition to the interaction Hamiltonian. This 
addition was necessary in order to satisfy the commuta- 
tion relation for time-dependent Bose-operators. If the 
interaction Hamiltonian has the normally ordered form 
|^| then it is not necessary to have deal with thermal 
fluctuations. This approach allows us to develop an al- 
gebra of time-dependent Bose-operators and to investi- 
gate the spectrum of quantum fluctuations of the quadra- 
ture components. The results of the quantum theory of 
self-action for USPs in the medium with the relaxation 



Kerr nonlinearity based on the normally ordered inter- 
action Hamiltonian and the developed algebra of time- 
dependent Bose-operators are presented below. 



II. THE QUANTUM EQUATION OF 
SELF-ACTION OF USPs 

For the monochromatic radiation the quantum equa- 
tion of self-action can be found, for example, in ||. Up 
to now, the proper quantum equation of the self-action 
for USPs with the account of relaxation behavior of the 
nonlinearity is absent in the literature. It is necessary 
to mention that the deduction of the quantum equation 
is based on the interaction Hamiltonian. However, in 
this case we obtain the time-evolution equation for the 
Bose-operators. For monochromatic radiation the con- 
version of time-evolution equation into space-evolution 
one is done using the t — >z/u replacement, where z is 
time variable and u is the speed of pulse in nonlinear 
media. If we deal with the propagation of pulse in non- 
linear media then the " impulse operator" of a pulse field 
should be used We begin with the analyse of the 
self-action of UPSs in non-incrtial nonlinear media. 



A. THE QUANTUM EQUATION OF 
SELF-ACTION IN THE NON-INERTIAL 
NONLINEAR MEDIA 

In nonlinear media with non-inertial behavior the self- 
action process is described using_the impulse operator 
(quantity of movement) Gi„t(z) 



G int {z)=f3h N[n 2 (t,z)}dt, 



(1) 
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where N is the operator of normal ordering, factor (3 is 
defined by the cubic nonlincarity of the medium In 
consequence, in the Hcisenberg representation the quan- 
tum space-evolution equation for the annihilation pho- 
tons Bose-operator in a given cross section z (A(t,z)) 
has the form 



ih 



dA(t, z) 
dz 



A(t,z),G int (z) 



(2) 



In agreement with (Q) the quantum equation of self- 
action for a light pulse follows from (0) 



dA(t, z) 
dz 



i(3A + (t 1 z)A 2 (t,z) = 0. 



(3) 



Eq.(g) is written in the moving coordinate system: z=z' 
and t=t' —z/u. It is important to note that in compar- 
ison with the so-called nonlinear Heisenberg equation, 
used in the quantum theory of optical solitons, in (J3[) the 
dissipation of light pulse in the nonlinearity is not taken 
into account. This approach corresponds to the first ap- 
proximation of dissipation theory. In fact, the traditional 
way to introduce the quantum equation of self-action is 
based on the interaction Hamiltonian. In this case one 
gets the time-evolution equation. The transition to the 
space-evolution, as already mentioned, is realized using 
replacement t — > z/u. This approach is enough reason- 
able in case the radiation is monochromatic. If we deal 
with the nonlinear propagation of a pulse, we use the im- 
pulse operator of a pulse field ([l ) which is connected with 
the evolution of field in space. Eq.(||) has the solution 



A(t,z) 



3(t) io(*), 



(4) 



where r )=f3z and, as usual, Ao(t) is the value of the op- 
erator at input of nonlinear media (Ao(t)=A(i, z=0)), 
•ho(t)=AQ (t)Ao(t) is the photon number "density" oper- 
ator. For its hcrmitian conjugated operator we find 



A+(t,z) = A+^e- 1 ^ ^. 



(5) 



In agreement with ([|) and (||) the operator h{t, z) ~ 
A + (t, z)A(t, z) does not change itself in nonlinear 
medium: 



hit, z) = h(t, z=0) = ho(t), 



(0) 



where z=0 corresponds to the input of the nonlinear 
medium. In fact, (|^) means that the photon statistic 
in media remains unchanged. The commutation rela- 
tion at the input (z=0) of nonlinearity [Ao(ti), ^o"(^2)] = 
5(^1 — ^2) , should be satisfied for any coordinate z in non- 
linear media: 



[A(t 1 ,z),A + (t 2 ,z)] = 5(h-t 2 ). 



(7) 



The solutions and (| 
commutation relation (m 



do not permit to verify the 
Besides, the analyse of the 



statistical characteristics of the pulse is accompanied by 
the necessity of the reduction to the normally ordered 
form of the expression e lln °^ . In this case, the solu- 
tions (|J) and (|^) are accompanied by the singularity of 
the function S(t) at t=0. The specified circumstances 
represent the main deficiency of the quantum theory of 
self-action of USPs in non-incrtial nonlinear media. 



B. THE QUANTUM EQUATION OF 
SELF-ACTION IN THE INERTIAL NONLINEAR 
MEDIA 

In the classical theory, the self-action process in inertial 
nonlinear media is described by the equation (in the first 
approximation of the dispersion theory) p] 



dB(t,z) ldB(t,z) 
dz u dt 



. k 



1— An{\B{t,z)\ 2 ) 
n 



xB(t,z)=0, (8) 



where: B(t, z)- the complex amplitude of a pulse, z- the 
distance in the nonlinear media, u- the group velocity, 
An(\B(t, z)\ 2 )=An((t, z)- the nonlinear addition to the 
coefficient of refraction. We consider that the last term 
of (|[) is caused by the high-frequency Kerr effect, and its 
evolution follows from the equation 



dAn(t, z) 
' dt 



1 



+ An(i,z) = -n 2 \B(t,z)\' 



(9) 



Here r r represents the relaxation time of the nonlinearity 
and 712- the nonlinear factor. We mention that in general 
the behaviour of the nonlinear addition differs from the 
one characterized by (||). However, if the carrying fre- 
quency of a pulse is far enough from the resonance and 
the pulse duration t p is greater that the relaxation time 
r r , then (H) is correct The solution of (||) looks like 



An(t, z) = (n a /2) f H(t - h)\B{t u z)\ 2 dh. 

J — OO 



(10) 



The function of nonlinear response H(t) is entered in 
such a way that in the limit r p 3>r r the nonlinear addi- 
tion becomes An(t, z)=(ri2/2)\B(t, z)\ 2 . In the moving 
system of coordinates (t'—t—z/u, z'—z), taking into ac- 
count dlQ), cq. (ph takes the form 



dB(t,z) 
dz 



poo 

ij* ff(ti)|B(t-ii 
Jo 



z)\ 2 dhBfaz), (11) 



where 7*=/co r W2no and the quotation-marks ( ' ) in new 
system of coordinates further will be lowered for simplic- 
ity. The transition to the quantum equation usually is 
carried out in the spectral representation. However, in 
the considered case it is more natural to use time repre- 
sentation. We make in ( |ll| ) replacement of the complex 
amplitudes with the operators, entered in the previous 
sections, 
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B{t,z)^iCA(t,z), B*{t,z) -> -iCA + (t,z). (12) 

The right part of the equation we have gotten in this way, 
will be written below in normally ordered form. In order 
to take into account the presence of the vacuum fluctua- 
tions, existing up to the moment of arrival of the pulse, 
we will replace in jl]] ) the bottom limit of integration by 
-oo. As a result we get (C=(?iwo/2^) 1/2 , /3=7*C 2 /2) 



i/3 / H{\t l \)A+{t-t 1 ,z)A{t-t 1 ,z) 



dA{t, z) 
dz 

xA{t,z)dti. (13) 

Eq. ( |l3| ) represents the correct quantum equation and can 
be obtained from space evolution equation for the opera- 
tor A(t,z) in interaction representation (^|). Taking into 
consideration the incrtial behaviour of the nonlincarity, 
the impulse operator of a pulse should be introduced as: 

/■OO ft 

G mt {z)=hf3 dt H(t-h) 



oo «/ — oo 



x'N[n(t,z)n{t 1 ,z)]dt 1 , (14) 

where H(t) is the nonlinear response of the medium (see 
©) (H{t)^0 at t>0 and H(t)=0 at t<0). We note that 
the integral expression in (UJ) at the moment of time t 
depends only on the previous ones. Therefore, in this 
case the causality principle relative to measured physical 
value is not broken. An important condition that must 
be satisfied by the impulse operator is 



[G int (z),n(t,z)] = [G mt (z),n (t)} =0, 



(15) 



which means that the photon number operator remains 
unchanged in medium (photon number is a constant of 
motion f[j]l)- The annihilation Bose-operator which agrees 
with the ( |l4| ) verifies the equation of self-action 



dA(t, z) 
dz 



ipq[n {t))A{t,z) = Q, 



where 



#o(t)]= / #(ti)[n (t-ii) + n (i + ti)]dti. 
Jo 



(16) 



(17) 



If H(t)=S(t) then @ is converted into @). Eq.(|) can 
be obtained in the same form if we consider that the re- 
sponse of the nonlincarity has relaxation behaviour H(t) 
in accordance with ( [To|) and the relaxation time i> is 
significantly less than the duration of a pulse t p . At 
the same time, as will be shown further, in this limited 
case also, the account of finite relaxation time plays the 
main role in formation of the nonclassical light. It is 
necessary to note, that by replacement in ( |l6| ) the oper- 
ators on complex values we obtain the equation, which 
not completely coincides with the classical equation of 
self-action in presence of the non-stationary nonlinear 
response. There is no second term, included in dlj]) as 



fio(t + ti). The presence of this term in quantum theory, 
which at the first sight is in contradiction with the causal- 
ity principle, is connected with the quantum description 
that even in absence of a pulse the vacuum fluctuations 
always are present. A response function was introduced 
already in the model for the Kerr effect by Blow et al. 
H . Although the need for an attendant noise source was 
anticipated by these authors, they did not indicate where 
it should be inserted. In |Q the quantum noise as ther- 
mal fluctuations was additive inserted in the interaction 
Hamiltonian but this procedure did not allow to develop 
the consistent quantum theory of self-action of USPs. Ac- 
cording to Jl4|), the operator n(t,z) remains unchanged 
in nonlinear media (see @). Solving the spatial evolu- 
tion equation (fl6l), the annihilation (creation) photons 
Bose-operators in nonlinear medium have the following 
form: 



A(t, z) = e^ q[AoW] A (t). 
A + (t,z) = A+(t) e 



-i-yq[h (t)] 



(18) 
(19) 



In dH) and (Q A (t)=A{t,0), y=f3z. The expression 
q[no(i)] (see (|l7[)) it is convenient to be written as: 



hit^hoit - h)dh, (h(t) = H(\t\). (20) 



If we consider h to be time-independent, then the ex- 
pressions (18) and ( fl9| ) give results for the monochro- 
matic field. In the case that the nonlinear response in 
(|l8|), ( |Tg| ) has the form h(t) = 8(t), the results for non- 
inertial nonlinear media can be obtained (sec @). To 
find the statistical characteristics of a pulse at the out- 
put of the nonlinear medium it is necessary to estimate 
the averages of the operators A(t,z), A + (t, z) and their 
combinations. They can be estimated, if the operator ex- 
pressions are given in the normally ordered form, when 
the creation photons Bose-operators are placed at the left 
of the annihilation photons Bose-operators. The use of 
the expressions (18), ( |l9| ) involves the development of a 
special mathematical device. 



III. THE ALGEBRA OF THE 
TIME-DEPENDENT BOSE-OPERATORS 

For the beginning, in order to simplify some expres- 
sions we introduce the operators: 

6{t) = nq[h (t)}, 6+(t) = -ijq[n (t)}, (21) 

where + (t)= — 0(t). Hence, the equations of self-action 
(Il8|), ([[9]) can be represented as: 

A(t,z)=e 6 ^A (t), (22) 
A+(t,z)=A+(t)e 6+ ( t \ (23) 

Taking into account (pO) , is it easy to remark that 
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[0(h),0(t 2 )} 
In consequence we have 



e o(ti) e o(t 2 ) 



e O(t 2 ) e O(t!) 

_<5(ti)+6(i 3 ) 



3 6(t 2 )+o(ti 



(24) 



(25) 



A. THE PERMUTATION OPERATOR 
RELATIONS 

The following operator permutation relations hold: 
i (*i)O(t 2 ) = [6(t 2 ) + ijh(t 2 - t^Aoih), 

A (h)6+(t 2 ) = [6+{t 2 ) - ijhih - h)]Mti), 

A+it^Oih) = [6(h) - i 7 ft(t a - ii)]iif(ti), (26) 
ij(ii)0+(i 2 ) = [0 + (t 2 ) + ijh(t 2 - 

Using the mathematical induction it is possible to demon- 
strate the validity of the formulae (m 6 AT): 

io(ii)O m (t 2 ) = [6{t 2 ] + z 7 /i(t 2 - t^rAoih), 
A (t 1 )6+ m (t 2 ) = [6+{t 2 ) - vyh(t 2 - *i)Pio(ti), 
i+(ti)0™(t a ) = [0(t 2 ) - i 7 M*2 - h^A+ih), (27) 

i+(t 1 )d+ m (t 2 ) = [d+(t 2 ) + i 7 h(t 2 - t^r A+ih). 

To simplify the operator algebra is useful to redefine 
G(t2-fi)=*7J»(*2-«i), G*(t 2 -ti)=-i7/i(i 2 -ii). (28) 

Decomposing e^'- 1 and e° + (*' in Taylor series we obtain 
the operator permutation relations which play an impor- 
tant role at the estimation of the statistical characteris- 
tics of a pulse. Hence, finally we have: 



e°^A+(t 2 )=A+(t 2 )e 



+ (t„\ P 6<,ti)+G(t 2 -t 1 ) 



(29) 



>? + ^A+(t 2 )=A+(t 2 )e 



+ /f„\ P d + (t 1 )+G'(t 2 -t 1 ) 



Using the permutation relations (|29j) it is possible to ver- 
ify the commutation relation (Q) for the operators A(t, z) 
and A + (t,z). 



B. THE NORMAL ORDERING THEOREM 

As pointed out previously, in the considered analyse 
another important question is represented by the reduc- 
tion to normally ordered form of the operators A(t, z) 
and A + (t, z). In (EG) we proceed to normalized time 
= ti/r r and then (|21| ) can be presented like: 



0(t) = 17 / h(6)h (t - 0r r )c 

-QO 

G{e)n (t-6r r )de, 



(30) 



where /i(0) = r r /i(0T r ) and G(0) = 17/1(0). For the re- 
duction of the expressions (|2l|) to the normally ordered 
form the following theorem can be used: 
Theorem: Bose-operator e '*' can be represented in 
the normally ordered form this way: 



Ncxp 



,G(9) 



n (t - 6T r )d6 \ . (31) 



The operators in the integral expression may be under- 
stood as the c-numbers. In ]l0| the similar theorem in 
the spectral representation was formulated and demon- 
strated and we mention only that the similar demonstra- 
tion can be done also in this case. In fact in jnj the 
integration limits are not defined and the theorem has 
not a obvious applicability. The demonstration of ( |3l| ) 
does not represent the central objective of this article so 
we formulate the theorem in the time-representation only. 
The average value of the is given by the formula 



oO{t) 



cxp 



n (t - 8r r )de \ . (32) 



C. THE AVERAGE VALUES OF e 6{t) AND e 6+w 

In most of the experimental situations 7 « 1 which 
allows one to decompose the integral expression in the 
( |32"| ) and to limit decomposition at terms having order 
7 2 . Using this approach we have: 



(e (5( * ) )=exp / G(0)n o (t - dT r )d9 



G 2 {9)n (t-6T r )d6 



(33) 



where fio(t)=(ho(t)}=\ao(t)\ 2 . It is convenient to enter 
in further analyse the envelope of a pulse p(t), so that 
o:o(t)=p(t)ao. If the initial pulse has the gaussian form 
then p(0)=l. For simplicity we denote: 



i(){t) = 7 / h{8)n (t- 8T r )d9, 
J —00 

1 f°° ~ 
=~1 2 h 2 (0)n o (t~OT r )dO. 



From 



we have 



a iil>(t)-ii(t) 



3 -ii/-(*)-M(t) 



(34) 
(35) 

(36) 
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The parameters ip(t) and p(t) are connected with the 
self-action effect and tp(t) represents nonlinear phase ad- 
dition. Then 



m = ^o h(9)p 2 (t-er r )de, 



fi(t) = no / h 2 (o) P 2 (t-0T r )d9, 

Jo 



(37) 
(38) 



where ipo = 270:9 = 2jfio and po = 7 2 "-o = 7"0o/2. A 
special interest is represented by the estimation of the 
average values of the Bose-opcrator combinations at co- 
herent initial states (see (pl|)). Taking into account (J2T 
we have: 



^OiH) 0(t 2 ) 0(t 2 ) Oi^) _0(t 1 )+0(t 2 ) _ 0(^,12) 



e O(*i) e O(* 2 ) 



= N 



?(*1,*2) 



In consequence we find 

N 
where: 

Q{t u t 2 ) = 6(h) + 6{t 2 ) 

G(h,t 2 ;6)no{6T r )d6, 
G(h,t 2 ;8) = ijh(h - 9T r ) + ijh(t 2 - 6r r ). 



(39) 



(40) 



(41) 
(42) 



Using the theorem of normal ordering for (41) we es- 
timate the averages of different combinations of Bosc- 
operators: 

/ e 6(ti)+6(t 2 )\ _ e i[V>(*i)+V'(*2)]-M(*i,*2)-if(ti,t2) ) 

(gO + (*l)+0(*2)\ _ e ll-4>(tl)+1p(t 2 )]-tl(tl,t2) + K(t u t 2 ) 



6(t 1 ) + 6+ (t 2 )\ _ i[Vl(tl)-V(«2)]-M(*l.*2)+if(il,t2) 



(43) 



(5+(ti)+6+(t 2 )\ 



-i[il){t 1 )+-4>(t 2 )]- t i(t 1 ,t 2 )-K(t 1 ,t 2 



where p(ti,t 2 )=p(ti)+p(t 2 ) and K(ti,t 2 ) represents the 
temporal correlator 

/oo 
h{t 1 ~eT r )h{t 2 ~eT r ) P 2 {9T r )de. (44) 
-oo 

In agreement with most of the experimental situations 
the a pp roximation r p 3>T r can be used. In this case in 
([}?]), (|3^) and (^J) we can suppose that in time the enve- 
lope of a pulse slowly change itself so it practically does 
not depend on the change of the integration variable. 
Therefore it is possible to eliminate it from the under in- 
tegral expression in the essential point 9r r = in ( ^7]|3^ ) 
and 9r r = t\+r/2 in (Q) which corresponds to the max- 
imal value of under integral expression h{6) = 1 (p7[ 3S) 
and h(t 1 -0T r )h(t 2 -dT r )=h(-T/2)h(T/2)=h 2 (T/2) (gj) 
consequently. Then 



ij(t) = ^oP z (t) / h(6)d6, 
Jo 

/>OC 

p(i) = p p 2 (t) / h 2 (6)d6, 
Jo 



(45) 
(46) 



K(t 1 ,t 2 ) = p p 2 (t 1 +T/2)— h(6)h(0 + T)d8. (47) 

T r J-oo 

Here T=t 2 ~t\. One should note that in our previous 
analyse we did not choose the relaxation function of non- 
linearity in a definite form. If the nonlinearity is of a 
Kerr type, the relaxation function should be introduced 

as 

#(f) = (l/7- r )e-*/ T - (t>0). (48) 
Then fr(6>)=e~ |e and for the integrals in we find 



~ f°°~ 1 

h(0)<W = 1, / h 2 (9)d9 = -, (49) 
o Jo 2 



g(r) = — [°° h{8)h(9+T)d6 = - fl+ti)h(- ) . 



IV. THE CORRELATION FUNCTION OF 
QUADRATURES 



As stated earlier (see (|6j)) in self-action process the 
photon statistics remains unchanged. Therefore, we are 
interested in analysing the quadrature components which 
are defined as: 



X{t,z) = A(t,z) + A + (t,z) /2, 
Y{t, z) = \A(t, z) - A + (t, zj\ /2i 



(51) 
(52) 



The averages of the operators A(t, z) and A + (t, z) at ini- 
tial coherent state of a pulse are 



(A(t,z)) =a (t)(e d ^), 
(A+(t,z)) = a* (t){e 6+ V). 



(53) 
(54) 



Taking into account ( p6[) and that ao(t)=\ao(t)\e l ' p ( t \ for 
average values of quadratures we obtain: 



(X(t,z)) = |oo(t)|e-"W cos*(t), 
(Y(t,z)) = |ao(*)|e-" (t) sin$(t). 



(55) 
(56) 



where $(i) = ip(t)+(p(t). Exponential term in (]5q ) and 
( f56| ) is caused by quantum effects - in the classical theory 
it is not present. From j55|)-(56) it is concluded that 
the changes of quadratures in time are connected with 
changes in pulse's envelope. We introduce correlation 
functions of quadrature components as 
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Rx{t, <+^)=2 P (*> «)Jf(t+T))+(A-(t+T, «pf(t, z)) 



2(X(i,z))(l(i+r,z)) 



(57) 



fl y (f, <+r)=- (F(t, Z )y( f+r ))+(y(t+ rj 2 )y(t, *)) 



-2(Y(t,z))(Y(t+r,z)) 



(58) 



To analyse the correlation functions of quadrature com- 
ponents, it is necessary to evaluate the correlators 
U(t 1 ,t 2 )={X(t 1 )X(t 2 )) and £ Y (ti,t2)=(Y(ti)Y(t2)). 
Using permutation relations (129) and (E3h, we obtain 



- ti) + iMii)IMt2)|e-^'' 2) 

e - A (*i>*»)c08[$(ti)+$(t 2 )-|-7fe(t 2 -ti)] 



A(*1,*2 



cos[$(ti)-$(t 2 ) 



(59) 



e~ A(tl '' 2) cos[$(i 1 )+$(i 2 )+7^(t 2 -t 1 )] 



= A(ti,t 2 ) 



cos[$(ti)-$(t 2 ) 



(60) 



where A(*i, i 2 )=/i(ii, t 2 )/i(i 2 — ii) (see @)). As a result 
for the correlation functions of quadratures we have: 



flx(t,t + r) 



i r 



S(t) 



R Y (t,t + T 



-^op(t)p(t+T)h(r) sin[$(i) + <P(t+r)} 
+i>oP(t)p(t+T)g(t, t) sin $(t) sin $(i+r) 

l r 



where 



+VoP(*V(i+ T )M T ) sin[$(t) + $(t+r)] 
+ip%p(t)p(t+T)g(t, T ) cos cos $(*+t) 

5 (t,T)=p 2 (t + r/2) 5 (r). 



(61) 



(62) 



(63) 



To obtain (|6l]) and ( p2| ) the 7 <C 1 and r r <C t p approxi- 
mations have been used. 



V. THE SPECTRUM OF QUANTUM 
FLUCTUATIONS OF QUADRATURE 
COMPONENTS 

Spectral densities of fluctuations of the quadratures are 
defined by the following expressions: 



/oo 
R x (t,t + r)e iur dT, 
-OO 

/OO 
R Y {t,t + T)e lulT dT. 
-OO 



(64) 
(65) 



Taking into account the weak change of the envelope dur- 
ing the relaxation time one obtains: 



1 - VoP 2 (*)sin2$(i) / h(r)e iuT dr 



+%p 4 (i)sin 2 $(t) / ff (r)e^ r dr , (66) 

1 

4 



S Y (u),t) = - 1 + VoP (*) sin2$(i) / h(T)e luJT dT 



-H/> 2 p 4 (i)cos 2 $(t) / g{T)e luJT d T . (67) 



The estimation of integrals in (|6(j) , (q7|) gives us 

2L(Sl), 

2 =4L 2 (r»), 



1 + (cJT r ) 2 

4 

[1 + (ur f ) 



/i(T)e MT dr 

O 

3 

g^e^dr 

where O = wr r . Hence 

S x (to,t) = - l-2V'(t)£(n)sin2$(i) 

+4^ 2 (t)L 2 (fi)sin 2 $(t) 
1 + 2^(f)i(0) sin2$(i) 
-f4^ 2 (i)i 2 (^)cos 2 $(t) 



(68) 
(69) 



5y(fi,t) 



1 r 



(70) 



(71) 



where ip{t)=2~/\a (t)\ 2 . From @ and (|l|) follows that 
the choice of the phase §(t) determines the level of quan- 
tum fluctuations lower and higher than the shot-noise 
level Sx(v)=S y (lj)=1/A, corresponding to the coherent 
state of the initial pulse. In conformity with the Heisen- 
berg relation the behaviour of the spectrum of the X- 
quadrature appears to be moved with a phase 7r/2 in 
comparison with y-quadrature. In case of an optimal 
phase of the initial pulse 



arctan 



i>(t)L(n ) 



(72) 



chosen for the frequency £lo=LUoT r , spectral densities 
1 



s x (j] ,i)=-[vi+?PW-^Wi(«o) 



s Y (n , t) =- yi+^{t)L 2 {n )+^{t)L{n ) 



(73) 
(74) 



Eqs . (|73|) , (|74|) indicate that when the nonlinear phase ad- 
dition ip(t) increases Sx(^o,t) monotonously decreases 
and S Y (p,o,t) monotonously increases. At any frequency 
Q we have 

s x (n,t) = s x (n ,t) + ^i>(t)[L(n) - L(n )] 

x{[L(fi)+L(fi )]V7(t)-[l+(L(r!)+ J L(r! )) J L(fi )V 2 (i)] 
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FIG. 1. Dispersion [Z] of the quadrature component of a 
pulse at time i=0 as a function of the maximum nonlinear 
phase ijj(0)=2'y\a 2 (0)| [X] and the reduced frequency fl=wr r 
[Y], at values of the phase of the pulse which are optimal for 

n=o. 



FIG. 2. Dispersion [Z] of the quadrature component of a 
pulse at time t=0 as a function of the maximum nonlinear 
phase V'(0)=27|a 2 (0)| [X] and the reduced frequency fl=W7y 
[Y], at values of the phase of the pulse which are optimal for 
0=1. 



:[l+^(t)L a (fi)]-Va}, 



(75) 



S Y (n, t) = S Y (Q ,t) + -^(f)[Z(fi) - L(f2 )] 

x{[L(r!)+L(f7o)]V'W+[l+(i(«)+i(«o))i(«o)^ 2 W] 
xfl+^^L 2 ^)]- 1 / 2 }, (76) 

The spectra of X-quadrature component, calculated 
by the formula (|75|), at t=0 (ip(0)=ipo) for the cases 
flo=0 (wo=0), Oq=1 (wo^^ 1 ) are presented in Figs. [l],|| 
respectively. On Fig. |l| one can see that for wq=Q spec- 
tral density of X-quadrature component is minimal on 
frequency lu=0 for any values of phase ipo- For ^ot^O 
(Fig. ^ and phases "0o>l the minimum of the fluctu- 
ation spectrum of X-quadraturc component lies at fre- 
quencies Q=l (u>=T r r 1 ), and for ?/>o<l the minimum lies 
near fl^O (w~0). 



VI. THE WIDTH OF THE SPECTRUM OF 
SQUEEZED QUADRATURE 

From Fig. [j] one can conclude that the frequency where 
spectral density of "-quadrature fluctuations" is lower 
than the shot-noise level, depends on the nonlinear phase 
addition %pQ. Width of the spectrum below the shot-noise 
level AQ=T r Aoj should be defined from 



1 



1 



S x (Afl,t) = -[-+S x (fl ,t) 
Accounting (J7|j7|) for Q a =0 from @ we have: 



(77) 



+ V(i)v / 1+W)-V' 2 (*)-1 = 0. (78) 

Eq.@ in L(Afi) (sec @|9|))has two solution of which 
only one is real. Solving the (ym for Afl we get: 



An 



^)Vl + ^ 2 (t)-V 2 W-l 



1/2 



(79) 



From ( |79| ) it follows that the change of Afl is connected 
with changes in pulse's envelope. The frequency band in 
which the spectral density of the quadrature fluctuations 
is lower than the shot-noise level depends on the non- 
linear phase shift tp(t). The corresponding dependence 
at t=0 for oj=0 is displayed in Fig. |j. It may be noted 
that at V'o^l width of the spectrum below the shot noise 
level is one and a half width of the spectral response of 
nonlincarity. 



VII. PHOTON NUMBER SPECTRAL "DENSITY" 
OF PULSES WITH SELF-PHASE MODULATION 

The spectral photon number operator is defined by 

h(ui, z) = a + (oj, z)a{u>, z), (80) 
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FIG. 3. Spectral band Afi=27yAtj of the quadrature of a 
pulse with suppressed quantum fluctuations as a function of 
the maximum nonlinear phase tpo. 



H(6) 



?v eXP 'l 2r 2 j 



(86) 



then K{t,t+r) = p p 2 (t+Tj2) exp {-t 2 /4t 2 }. Using the 
expression K(t, t + r) in (|83|) in paraxial approximation 
we have 



n(f2, z) 



n 



1/2 



(l + 7V'o^ 2 /4) +4^ 2 

f! 2 (1 + 7^74) 



x exp 



(l+ 7 ^ 2 /4) 2 +4^ 



(87) 



where v—T p /r r . From ( |S7| ) one can conclude that the 
spectral density depends on the nonlinear phase addition 
and on the relation between the pulse duration and the 
relaxation time of the nonlincarity. 



VIII. THE CORRELATION FUNCTION OF 
SPECTRAL COMPONENTS OF SPM-USPs 



where: 



a(u>, z) 
a + (u>, z) 



2ir t„ 



2tt t. 



A(t, z)e luJt dt, 
A+{t lZ )e- wt dt. 



(81) 
(82) 



p J — oo 



Taking into account ( |S1,[82"| ) and (|22| , |23| ) for photon num- 
ber spectral "density" (8(3) we find: 



n(u>, z) 

where 
I = 



2^r 2 



OO /» OO 



oo J — oo 



x e M*»-*0 (ftl dt 2 s5°.-L.|I 
2tt t 2 



(83) 



p(t) exp {i [ip p 2 {t) + ip(t) + cut] }dt. (84) 



The last expression in ( 83 ) is written without the account 
of p:(t\,t 2 ) and K(t\ , t 2 j[see (fl3|)). If the initial pulse has 
gaussian form p{t)= exp {— i 2 /2r 2 } then, using paraxial 
approximation (p 2 (t) = 1— t 2 /r 2 ) JllJ in ( [S4"|) for spectral 
density (B3h we find 



W>class (^j z) 



"0 



Vl+4^ 



exp 



O 2 



1 + 4^ 



(85) 



where Q=lut p . From (|85J) it follows that the spectral 
density of a pulse with self-phase modulation (SPM-USP) 
decreases when nonlinear phase addition increases. To 
calculate ( |35"| ) in ( |S3| ) the terms /u(ii,t 2 ) and K(ti,t 2 ) 
have not been taken into account. In consequence, the 
spectral density (|8^) does not depend on relaxation time. 
If we take the relaxation function of nonlinearity as [O] 



We introduce the correlation function of different spec- 
tral components in the following symmetric form: 

R(wi, u> 2 , z)=- {h{u}\, z)n(uj 2 , z)) + (h(uj 2 ,z)n(uj 1 , z)) 



-2(n(uii, z))(h(ui 2 , z)) 



(88) 



Leaving out the preliminary accounts for (88) we obtain 



fl(«i,W2, z) = hSfo-wi) - ^ Im{I*I 3 + I 2 i;}, (89) 



h = 



OO /"OO 



p(t 1 ,t 2 )exp{i^ 2 (<i)- ^ 2 (t 2 )}} 



— oo J — oo 



X 



expz[uiti — uj2t2]dtidt2, (90a) 



OO /*OC 



h = I / p(h,h)e^{i[^(h) + ^ 2 (t 2 )}} 

-oo J — oo 

x exp{i[u;iti +^ 2 ^]}^i^2, (90b) 



OO /"OO 



Js = / / g(ti,i 2 )exp{i[^(ii)+V (* 2 )]} 

J — oo J — oo 

x exp{i[ujiti + Lu 2 t 2 ]}dtidt 2 . (90c) 

?(fi,ta)=p(*i)/>(*a), q{tiM)=p{h,t 2 )h{t2-tx). If the 
initial pulse has gaussian form and the relaxation func- 
tion has the form (|8^) then in paraxial approximation 
for correlation function (p9[) we find 



i?(0 1 ,Q 2 ,z)=/ 1 5(Q 1 -Q 2 ) 



00 



2cW/3 



/m{r(Oi,0 2 )}, (91) 



h = [n das ,(» 1 )n c ; Q , s (» 2 )] 1 / 2 exp^^ 1 2 L (92) 

L 1 + 4?/>n J 
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In fl9_y) are entered the following designations: 



a=[l + 4^ 2 ] 



1/2 



,2\2„/,21 1 / 2 



/3 = [1 + 2v 2 - 4t/> 2 ) 2 + 16(1 + u 2 f^l 



r(fii,0 2 ) = exp{G + iS} + cxp{E + iF}, 



Jmr(Oi,f2 2 ) = e G sinS , + e is sin.F, 



(93) 



(94) 



(95) 



(96) 



G = _Ol+j4 cose _^l cosS-^S cos(S-e) 
2a 2g 2/3 V ^ 

Vcos£ 2 -^cos(E + 0, 



/3 

S1 2 + S1 



2f?/3 



(97a) 



sin e H sin E 

2a 2g 2/3 



sin(E-0 



sin C+^i sin (E+0-(c+f ), (97b) 



2g/3 



n? + n% n 2 

E = 1 - 2 cos e cos E 

2a 2g 



^ cos (E-0 
2/3 



2,, 4 



F = 



rJ^v 2 cos^-^i^ cos (£+£), 
(3 2g(3 

n 2 + n 2 n 2 g n 2 , , 

sin e H 3- sin E — ^3- sin (E-£) 



(97c) 



2a 2g 2/3 

^ I/ 2 B i n f + ^ sin (S+0-(e+|) ) (97d) 
/3 2g/3 2 



g =[(! + , 2 ) 2 +^ 2 ] 1/2 , 
e = — arctan {2-00 }, 

4(1 + v 2 )^ 



£ = arctan 
E = arctan 



i + 2v 2 - m / ' 

l + i/ 2 



(98) 
(99) 

(100) 
(101) 



We define the spectral correlation function of the photons 
with frequency Q in the spectral band AQ 

Jn-An/2 Jn-An/2 

-no- (102) 

As a consequence, the conclusion which one can make is 
that for R(U,z)<0 take place the photon antibunching 
and for R(£l, z)>0 the photon bunching. 

The graphic dependence of the spectral correlation 
function fll02]) on ip at 0=0.04, Arj=2.5-10~ 3 and ^=10 
is displayed in Fig.||, whence it follows that the photon 




FIG. 4. Spectral correlation function Rasi(Q,z) [R] as a 
function of the maximum nonlinear phase ipo measured at 
the frequency fi=0.04 in spectral band Af2=2.5-10 -3 . 



bunching or antibunching can take place, and for phases 
ipo>l it becomes signi fican t. At frequency £1=0, the cor- 
relation function (see (102)) has the following simplified 
form: 



i?Afi(0,z) 



i/'o 



2&V/3 



(Aft) 2 sin 



-I 



(103) 



I 



and its dependence on tf>Q at Af2=0.75 is shown in Fig. 

whence it follows that the minimum of the spectral 
correlation function lies near ipoml. In this case the pho- 
ton antibunching takes place for all phases ?po>0 and it is 
maximal for tpQfal. It may be mentioned that the greater 
is the spectral band of measurement the stronger is the 
photon bunching or antibunching. 



DISCUSSION AND CONCLUSIONS 

The results presented in the present paper can be used 
for the correct interpretation of the results of experiments 
iSHIlllDl ' ^ n wmc h the laser pulses with the duration of 
the order 100 ps and quartz optical fibres were used and 
the maximal meaning of nonlinear phase shift ipQ was 
greater than 1. Certainly, in the measurement of the 
quadrature spectrum the suppression of quantum fluctu- 
ations of a pulse will be smoothed out (see (|75|)). This 
time over which the "smoothing out" occurs in the case 
of balanced homodyne detection JlJ] is determined by 
the duration of the heterodyne pulse. 

The developed theory enables the choice of the opti- 
mal strategy at producing and registration of ultrashort 
pulses in a squeezed states. The measurement of quan- 
tum fluctuation of short pulses take place at high frequen- 
cies of the order of several tens MHz in order to avoid 
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FIG. 5. Spectral correlation function Raq(£1, z) [R] as a 
function of the maximum nonlinear phase measured at 
the frequency = in spectral band Af2 = 0.75. 



any effects due to technical fluctuation concentrated at 
low frequencies. However, in this area the suppression of 
quantum fluctuations is greatest. The presented results 
show that by adjusting the phase of the signal pulse (or 
the phase of a heterodyne pulse) , maximal suppression of 
the quantum fluctuations can be realized at the spectral 
component of interest for us. This spectral component 
of interest can lie on the wing of the spectral response 
of nonlinearity (Fig. ||). This means, that for obtaining 
squeezed-light pulses the nonlinear media with a longer 
relaxation time and consequently with the greater non- 
linearity can be used [^|. 

Our results suggest that in the spectral measurements 
the photon antibunching can be observed. Usually, in 
the experiments spectral devices with confined spectral 
bands are used, thus limiting the amplitude of vacuum 
fluctuations which participate in the measurements. The 
final results show that the spectral correlation function 
depends on the nonlinear phase addition, the relation be- 
tween the pulse duration and the relaxation time of the 
nonlinearity and also on the spectral band of the mea- 
surement. In consequence the choice of the width of the 
spectral band of measurement can represent an effective 
method of control of the photon bunching or antibunch- 
ing. The obtained results indicate that the photon anti- 
bunching can be observed at any value of nonlinear phase 
addition in the low frequency measurements. At high fre- 
quency measurements the photon bunching or antibunch- 
ing strongly depends on the nonlinear phase additions. 

We note that the approach developed in the present ar- 
ticle can be used to analyse the formation of polarization- 
squeezed light in media with a cubic nonlinearity. This 
will be treated in a future publication. 
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